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Abstract. Propagation of Alfven waves in the solar plasma has been a topic of scientific interest 



o. 

for a long time. We have derived a dispersion relation lo 4 + lo 2 [(u 2 + rj 2 )k 4 — v\k 2 ] + iio(u + 
r])v 2 A k 4: + (vrjv\k e + i/ 2 rj 2 k & ) = for shear Alfven waves in a viscous and diffusive plasma. 
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The MHD equations for viscous and diffusive plasma are 



^ ■ p— — h p(v .V) v = — (Vx B)x B +pvV 2l v Momentum equation (1) 

V^H' OX P 

O 1 dB -* -» -» 

' -—— = V x ( v x B) + r]V 2 B Induction equation (2) 

_ 

V. 5= Magnetic flux conservation (3) 

where v is the velocity, B the magnetic field and p, p, rj, v are, respectively, the mass density, 
I/"") ' magnetic permeability, magnetic diffusivity and the coefficient of viscosity. Let us consider small 



perturbations from the equilibrium [1]: 

p = po + pi v = vt B=Bo + Bi 



O. 

t"-- , and linearize the equations (1) through (3) by neglecting squares and products of the small 

quantities (denoted by subscript 1). Here, the quantities with subscript represent their values 
in equilibrium. Notice that the equilibrium velocity Do is zero [1]. After linearization, we have 
the corresponding relations as the following: 



Po^- = -(VxBi)xB +p vV 2 v 1 (4) 

-^l = Vx(«ixBo)+']V 2 Bi (5) 

V.5i=0 (6) 

Here, the magnetic field Bo is taken uniform as well as time independent. On differentiating 
equation (4) with respect to time and then substituting the expressions for dv\/dt and dBi/dt 
from equations (4) and (5), we get 



(V x [V x («! xB )]| xB + -9_V x (V 2 Bi x Bo) 

>- > 1 1, On 



d 2 vi _ 2 

dt 2 L > ppo 

+— V 2 [(Vx Bi)x Bo] + v 2 V A v± (7) 
Wo 
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where Bq is the unit vector along Bo and va (= Bo/^fipo) is the Alfven velocity. Let us seek a 
plane-wave solution of the form 



Vi = v e 



i(k-r-uit) 



Bi=B e 



i(k-r-uit) 



(8) 



where k is the wave vector and lo the frequency. The effect of the plane-wave assumption is 
simply to replace d/dt by — iu; and V by ik- Under the plane-wave assumption, equation (7) 
reduces to 

u 2 V!= v\{ k x[k x(«i xB )]\ x B + — (fc . Bo) Bi 

ii//c 2 



+ [(fc x Bi)x Bo] ~ v k «i 

Wo 



For plane-wave assumption, equation (6) gives 

I • 5i= 



(9) 



(10) 



It shows that the magnetic field perturbation is normal to the propagation vector k ■ Equation 
(9) can be simplified as 

u 2 v 1= v\[(k -Bo){k -Bo) «i -(k .B )(k . v 1 )B + k [(k . «i) - {Bo- vi)(k -B )] } 

+—(k ■ Bo) Bi +— \(Bi (is ■ Bo)- k (Bo ■ Bij\ - v 2 k A v 1 (11) 
Wo Wo L J 

Let the propagation vector k makes an angle 6 with the equilibrium magnetic field Bo- Then 
equation (11) reduces to 

J 2 ~v 1= v A $ ] k 2 cos 2 9 vi -k cos 6(k . : v 1 )B + k (k ■ ~v\) - k cos 8(B . vi) } 



ir)k 5 Bocos 9 B\ wk 



+ 



k Bo cos 9 Bi - k (Bo ■ Bi) 



- v 2 k A ~v\ 



Wo Wo 
Multiplying equation (12) by Bo, we get 

Wo 

The imaginary part disappears when we consider B\ orthogonal to Bo- Now, we get 

(uj 2 + u 2 k A )(B . «i) =0 
Since (a; 2 + z/ 2 /c 4 ) is not zero, we have 

Bo- «i= 



(12) 



(13) 



(14) 



It shows that the perturbation of velocity is normal to the ambient magnetic field. The discussion 
so far leads to the geometry as shown in Figure 1. 
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Figure 1: Bo, Bi and «i are orthogonal to each other, k lies in the (Bo, «i)-plane and fc makes and 
angle 9 with Bo- 

Now, multiplying equation (12) by k and using equations (10) and (14), we get 
to 2 (vi . k) = v 2 A k 2 (v 1 . k) - v 2 k A (v 1 . k) 



u) 2 - v\k 2 + ^/c 4 (k . vi) = 



2 r A 



(15) 



Thus, we have either 



or 



k • v : 



oj 2 - v\k 2 + v 2 k 4 = 



(16) 



(17) 



The relation (16) leads to the shear Alfven waves whereas the relation (17) leads to the com- 
pressional Alfven waves. For shear Alfven waves, the perturbation is incompressible so that we 
have V. v i= 0. Under the plane-wave assumption, this relation gives the equation (16). 

Since k ■ «i= 0, we have 9 = 0. Using equation (16) (with 6 = 0) in (12) (for Figure 1), we 

get 



uj a vi= v\k 2 vi 4 



i(iy + r])k 3 B 
Wo 



Bi -v 2 k A v 1 



For the plane wave solution, equation (5) gives 

(rjk 2 — iuj) Bi= iBnk v\ 
From equations (18) and (19), we get a dispersion relation 

uj a + u?[(y 2 + 7] 2 )k 4 - v\k 2 \ + iw(u + V )v 2 A k 4 + [yr\v 2 A k % + v 2 V 2 k 8 ) 

This dispersion relation incorporates the viscosity as well as diffusivity. 
For the plane wave solution, equation (4) gives 

iB k 



(18) 



(19) 



(20) 



(yk 2 



10J) V\- 



Bi 



(21) 
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From equations (18) and (21), we get a dispersion relation 



2 r 2 
LO = k 



v\ — iuj{v + rj) + isr/k 4 



This dispersion relation is the same as obtained by Pekiinlu et al. [2] and can also be obtained 
by using the equations (21) and (19). 

We have thus obtained a new dispersion relation (20). 
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